DECOMPOSITION OF POLYNOMIALS AND 
APPROXIMATE ROOTS 



ARNAUD BODIN 

Abstract. We state a kind of Euclidian division theorem: given 
a polynomial P(x) and a divisor d of the degree of P, there exist 
polynomials h(x), Q(x), R(x) such that P{x) = ho Q(x) + R(x), 
with degh = d. Under some conditions h, Q, R are unique, and Q 
is the approximate d-root of P. Moreover we give an algorithm to 
compute such a decomposition. We apply these results to decide 
whether a polynomial in one or several variables is decomposable 
or not. 



1. Introduction 

Let A be an integral domain (i.e. a unitary commutative ring without 
zero divisors). Our main result is: 

Theorem 1. Let P G A[x] be a monic polynomial. Let d ^ 2 such that 
d is a divisor ofdegP and d is invertible in A. There exist h,Q,RE 
A[x] such that 

P{x) = hoQ( x )+R(x) 
with the conditions that 

(i) h, Q are monic; 

(ii) deg h = d, coeS(h, x^ 1 ) = 0, deg R < deg P - 
(hi) R(x) = Yl,i r i xl w ^h (deg Q\i =J> = 0). 

Moreover such h, Q, R are unique. 

The previous theorem has a formulation similar to the Euclidian 
division; but here Q is not given (only its degree is fixed); there is 
a natural Q (that we will compute, see Corollary 2) associated to P 
and d. Notice also that the decomposition P(x) = ho Q(x) + R(x) is 
not the Q-adic decomposition, since the coefficients before the powers 
Q l {x) belong to A and not to A[x]. 
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Example. Let P(x) = x 6 + 6x 5 + 6x + 1 £ Q[x]. If d = 6 we find 
the following decomposition P(x) —ho Q(x) + R(x) with h(x) = 
x 6 -15x 4 + 40z 3 -45x 2 + 30a;-10, Q(x) = x + 1 and = 0. If d = 3 
we have h(x) = x 3 + 65, = x 2 + 2x — 4 and = 40x 3 — 90x. 
If d = 2 we get h(x) = x 2 - ™, Q(x) = x 3 + 3ie 2 - \x + ^ and 
= -^a; 2 + 2|5 x. 

Theorem 1 will be of special interest when then ring A is itself a 
polynomial ring. For instance at the end of the paper we give an 
example of a decomposition of a polynomial in two variables P(x, y) £ 
A[x] for A = K[y\. 

The polynomial Q that appears in the decomposition has already 
been introduced in a rather different context. We denote by v^P the 
approximate d-root of P. It is the polynomial such that (\fP) d approx- 
imate P in a best way that is to say P — (\fP) d has smallest possible 
degree. The precise definition will be given in section 2, but we already 
notice the following: 

Corollary 2. 

q = Vp 

We apply these results to another situation. Let A = K be a field 
and d ^ 2. P £ K[x] is said to be d- decomposable in K[x] if there exist 
h,Q £ if [x], with deg h — d such that 

P(x) = hoQ(x). 

Corollary 3. Let A = K be a field. Suppose that char K does not 
divide d. P is d- decomposable in K[x] if and only if R = in the 
decomposition of Theorem 1. 

In particular, if P is d-decomposable, then P = h o Q with Q = \[P. 

After the first version of this paper, M. Ayad and G. Cheze com- 
municated us some references so that we can picture a part of history 
of the subject. Approximate roots appeared (for d = 2) in some work 
of E.D. Rainville [9] to find polynomial solutions of some Riccati type 
differential equations. An approximate root was seen as the polynomial 
part of the expansion of P(x)d into decreasing powers of x. The use of 
approximate roots culminated with S.S. Abhyankar and T.T. Moh who 
proved the so-called Abhyankar-Moh-Suzuki theorem in [1] and [2]. For 
the latest subject we refer the reader to an excellent expository article 
of P. Popescu-Pampu [8]. On the other hand Ritt's decompostion the- 
orems (see [10] for example) have led to several practical algorithms to 
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decompose polynomials in one variable into the form P(x) = ho Q(x): 
for example D. Kozen and S. Landau in [6] give an algorithm (refined 
in [5]) that computes a decomposition in polynomial time. Unification 
of both subjects starts with P.R. Lazov and A.F. Beardon ([7], [3]) for 
polynomials in one variable over complex numbers: they notice that 
the polynomial Q is in fact the approximate d-root of P. 

We define approximate roots in section 2 and prove uniqueness of the 
decomposition of Theorem 1. Then in section 3 we prove the existence 
of such decomposition and give an algorithm to compute it. Finally 
in section 4 we apply these results to decomposable polynomials in 
one variable and in section 5 to decomposable polynomials in several 
variables. 

2. Approximate roots and proof of the uniqueness 

The approximate roots of a polynomial are defined by the following 
property, [1], [8, Proposition 3.1]. 

Proposition 4. Let P G A[x] a monic polynomial and d ^ 2 such that 
d is a divisor of deg P and d is invertible in A. There exists a unique 
monic polynomial Q G A[x\ such that: 

deg(P-Q d ) <degP-^§^. 

o 

We call Q the approximate d-root of P and denote it by \fP. 
Let us recall the proof from [8]. 

Proof. Write P(x) = x n + aix n ~ l + a 2 x n ~ 2 + . . . + a n and we search an 

equation for Q (x) = x^ +b±x ! 3~ 1 + b 2 x 1 s~ 2 -\ h&«. We want deg(P — 

Q d ) < deg P— ^§^, that is to say the coefficients of x n , x n_1 , . . . , x n ~^ 
in P — Q d equal zero. By expanding Q d we get the following system of 
equations: 



{S) I 



a 2 = db 2 + Qbj 

a k = db k + c h-i k -i b i ■ ■ ■ frfc-i , 1 < k < 

ii+2i 2 H h(k— l)ik-i=k 



d 



where the coefficients Cj 1 ...j fc _ 1 are the multinomial coefficients defined 
by the following formula: 

d \ d\ 

C U---«fc-l 



ii,..., ik-ij i\\ - ■■ ik-iKd ~h h-i) } - ' 
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The system (S) being a triangular system, we can inductively compute 

the bi for % — 1, 2, . . . , ^: b x = b 2 = " 2 ^\ . . . Hence the system 
(iS) admits one and only one solution b\, b 2 , ■ ■ ■ , b*. 

Notice that we need d to be invertible in A to compute bi. Moreover 
bi depends only on the first coefficients ai, a 2 , ■ ■ ■ , a™ of P. □ 

Proposition 4 enables us to prove Corollary 2: by condition (ii) of 
Theorem 1 we know that deg(P - Q d ) < degP - so that Q is 

the approximate a?-root of P. Another way to compute \fP is to use 
iterations of Tschirnhausen transformation, see [1] or [8, Proposition 
6.3]. We end this section by proving uniqueness of the decomposition 
of Theorem 1. 

Proof. Q is the approximate d-ioot of P so is unique (see Proposition 
4 above). In order to prove the uniqueness of h and R, we argue by 
contradiction. Suppose hoQ + R — h'oQ + R' with R ^ R'; set rix 1 
to be the highest monomial of R(x) — R'(x). From one hand x l is a 
monomial of R or R', hence degQ \ i by condition (iii) of Theorem 1. 
From the equality (h' — h) o Q = R — R' we deduce that % = deg(R — R') 
is a multiple of deg Q ; that yields a contradiction. Therefore R = R', 
hence h = ti . □ 

3. Algorithm and proof of the existence 
Here is an algorithm to compute the decomposition of Theorem 1. 

Algorithm 5. 

• Input. P E A[x], d\ deg P. 

• Output. h,Q,R G A[x] such that P = ho Q + R. 

• 1st step. Compute Q = \f~P by solving the triangular system 
(S) of Proposition 4- Set hi(x) = x d , Ri(x) = 0. 

• 2nd step. Compute P 2 = P — Q d = P — hi (Q) — Ri . Look for 
its highest monomial OiX 1 . If deg Q\i then set h 2 (x) = hi(x) + 

i 

aiX de sQ , R 2 = Ri. If degQ \ i then R 2 (x) = Ri(x) + 0{X % , 
h 2 = hi. 

• 3thd step. Set F 3 = P — h 2 (Q) — R 2 , look for its highest 
monomial OiX 1 ,. . . 

• . . . 

• Final step. P n = P — h n -i(Q) — R n -i = yields the decom- 
position P = h o Q + R with h = h n _i and R = R n -i- 

The algorithm terminates because the degree of the Pj decreases at 
each step. It yields a decomposition P = h o Q + R that verifies all 
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the conditions of Theorem 1: in the second step of the algorithm, and 
due to Proposition 4 we know that i < degP — That implies 

coeS{h 2 ,x d ~ l ) = and degR 2 < degP - Therefore at the end 

coeff(/i, x d ~ v ) = 0. Of course the algorithm proves the existence of the 
decomposition in Theorem 1. 

4. Decomposable polynomials in one variable 

Let K be a field and d ^ 2. P G K[x] is said to be d-decomposable 
in K[x\ if there exist h,Q E K[x], with deg/i = d such that 

P(x) = /ioQ(x). 

We refer to [4] for references and recent results on decomposable poly- 
nomials in one and several variables. 

Proposition 6. Let A = K be a field whose characteristic does not 
divide d. A monic polynomial P is d-decomposable in K[x] if and only 
if R = in the decomposition P = h o Q + R. 

In view of Algorithm 5 we also get an algorithm to decide whether 
a polynomial is decomposable or not and in the positive case give its 
decomposition. 

Proof. If R — then P is ^-decomposable. Conversly if P is d- 
decomposable, then there exist h,Q G K[x] such that P = h(Q). As P 
is monic we can suppose h, Q monic. Moreover, up to a linear change 
of coordinates x \— > x + a, we can suppose that coeff (h, x^ 1 ) = 0. 
Therefore P = h(Q) is a decomposition that verifies the conditions of 
Theorem 1. □ 

Remark. Let P(x) = x n + a\x n ~ l H — • + a n , we first consider a±, . . . ,a n 
as indeterminates (i.e. P is seen as an element of K(ai, . . . , a n )[x]). The 

coefficients of h(x), Q(x) and R(x) = r x k + r 1 x k ~ 1 -\ \-r k (computed 

by Proposition 4, the system (S) and Algorithm 5) are polynomials in 
the dj, in particular rj = rj(ai, . . . , a n ) G K[ai, . . . , a n ], i — 0, . . . , k. 

Now we consider * G K as specializations of a i: . . . , a n and 

denote by P* the specialization of P at a*, . . . , a*. Then, by Proposi- 
tion 6, P* is (^-decomposable in K[x\ if and only if r^a*, . . . , a*) = 
for alH = 0, . . . , k. It expresses the set of (^-decomposable monic poly- 
nomials of degree n as an affine algebraic variety. We give explicit 
equations in the following example. 

Example. Let K be a field of characteristic different from 2. Let P(x) = 
x 6 +aiX 5 +a2x' k +a 3 x 3 +a4X 2 +a 5 x+aQ be a monic polynomial of degree 6 
in K[x] (the ai <E K being indeterminates). Let d = 2. We first look for 
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the approximate 2-root of P(x). \J P(x) = Q(x) 
In view of the triangular system (<S) we get 



x 3 + bix 2 + b 2 x + b 3 . 




a 3 - 2fri&2 
2 



Once we have computed Q, we get h(x) = x 2 + a 6 — b\. Therefore 

R(x) = (a 4 — 2bib 3 — b\)x 2 + (a 5 — 2b 2 b 3 )x. 

Now P(x) is 2-decomposable in K[x] if and only if R(x) = in if [x] 
that is to say if and only if (oi, . . . , a 6 ) satifies the polynomial system 
of equations in ai, . . . , a 5 : 



5. Decomposable polynomials in several variables 

Again if is a field and d ^ 2. Set n ^ 2. P G if [xi, • • • , x n ] is said 
to be d-decomposable in if [xi, . . . , x n ] if there exist Q G if [xi, . . . , x n ], 
and /1 G if [t] with deg h = d, such that 



Proposition 7. Let A = if [x 2 , . . . , x„] 7 P G A[x\\ = if [xi, . . . , x n ] 
monic in x\. Fix d that divides deg Xi P, such that char if does not 
divide d. P is d-decomposable in if [ ] if and only if the de- 

composition P = hoQ + Rof Theorem 1 in A[x\} verifies R = and 
h G K[t] (instead of h G if [t, x 2 , . . . , i n ]j. 

Proof. If P admits a decomposition as in Theorem 1 with R = and 
h G if [t] then P = h o Q is d-decomposable. 

Conversly if P is ^-decomposable in iff 1 then P — h o 

Q with /i G if[i], <5 £ if [xi, . . . , x n ]. As P is monic in xi we may 
suppose that h is monic and Q is monic in x±. We can also suppose 
coeff (h, = 0. Therefore h, Q and R := verify the conditions 
of Theorem 1 in A[x]. As such a decomposition is unique, it ends the 
proof. □ 

Example. Set A = K[y] and let P(x) = x 6 + aix 5 + (X2X 4 + a^x 3 + 
a^x 2 + (Z5X + a6 be a monic polynomial of degree 6 in A[x] = if [x, y], 
with coefficients = a,i(y) e A = K[y\. In the example of section 4 
we have computed the decomposition P = /io(5 + Pfor(i = 2 and set 



b, = f , b 2 = 22=£ 63 = ° 3 ~f 162 ■ We found h(t) = t 2 +a 6 -b 2 G A[t] 



a 4 — 2&X&3 — 62 
a 5 - 2b 2 b 3 = 0. 



P(xi, . . . ,x n ) = hoQ(x 1 ,...,x n ). 
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and R(x) = (04 — 26i6 3 — b\)x 2 + (05 — 26 2 6 3 )x G A[x] . By Proposition 7 
above, we get that P is 2-decomposable in K[x,y] if and only 

a 6 - 6§ G if, 

04-26163-62 = inX[y], 
a 5 - 26 2 6 3 = ini^y]. 

Each line yields a system of polynomial equations in the coefficients 
G K of P(x,y) = 'Y^,a i jX' l yi G y]. In particular the set of 
2-decomposable monic polynomials of degree 6 in K[x, y] is an affinc 
algebraic variety. 
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